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ABSTRACT 

Water  tunnel  tests  are  described  of  a  NACA  16-309 
section  equipped  with  a  25%  flap-chord  ratio,  simple,  sealed 
flap.  Lift  and  pitching  moment  coefficients  are  given  for  a 
range  of  angles  of  attack,  flap  angles  and  cavitation  numbers. 
Drag  data  are  not  considered  reliable  and  only  a  few  values 
are  presented.  The  results  are  also  given  of  limited  tests 
for  the  dependency  of  force,  moment  and  cavitation  inception 
data  on  Reynolds  number  over  the  range  from  1.25  to  3.98  x  106. 
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SOMMAIRE 

On  decrit  les  essais  qu'on  a  fait  subir  en  galerie 
hydraulique  a  une  section  du  profil  aerodynamique  NACA  16-309 
munie  d’un  volet  hypersustentateur  etanche  simple  ayant  un 
rapport  volet/corde  de  profil  de  25%.  On  presente  les  co¬ 
efficients  de  portance  et  du  moment  de  tangage  pour  une  gamme 
d'angles  d'attaque,  d'angles  d ' inclinaison  du  volet  et  de 
nombres  de  cavitation.  On  estime  que  les  donnees  relatives  a 
la  trainee  ne  sont  pas  fiables  et  on  n'en  donne  que  quelques 
valeurs.  On  presente  egalement  les  resultats  d ' essais  limites 
ayant  trait  a  la  dependance  de  la  force  de  sustentation,  du 
moment  et  de  la  cavitation  sur  une  gamme  de  nombres  Reynolds 
allant  de  1.25  a  3.98  x  106. 
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NOMENCLATURE 


section  drag  coefficient 

C  ^  section  lift  coefficient 

rate  of  change  of  lift  coefficient  with  angle  of 
attack  (degrees-1) 

C. ,  rate  of  change  of  lift  coefficient  with  flap  angle 

(degrees-  1  ) 

C  section  quarter  chord  pitching  moment  coefficient 

m 

rate  of  change  of  pitching  moment  with  flap  angle 
(degrees-  1 ) 

Re  Reynolds  number 

a  angle  of  attack  (degrees) 

6  flap  angle  (degrees) 

a  cavitation  number 

f  flap  chord  ratio 

x  abscissa  of  point  on  section  surface 

y  ordinate  of  point  on  section  upper  surface 

y^  ordinate  of  point  on  section  lower  surface 
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1.  INTRODUCTION 


The  NACA  16  Series  of  airfoil  sections  was  designed 
originally  for  high  speed  airfoil  applications,  where  a  uniform 
distribution  of  pressure  is  desired  on  the  upper  surface  at  de¬ 
sign  lift  coefficient.  This  condition  also  gives  favourable 
cavitation  characteristics,  a  factor  of  great  importance  in 
hydrofoil  section  design  since  cavitation  sets  practical  limits 
to  the  maximum  speed,  take-off  speed  and  degree  of  control  which 
can  be  obtained.  Consequently,  the  section  has  found  consider¬ 
able  use  in  fully-submerged  hydrofoils,  both  with  incidence 
control  and  with  a  simple,  sealed  flap  for  control  purposes. 

Results  were  published  in  1948  of  wind  tunnel  tests 
on  a  series  of  Type  16  sections1,  covering  a  fairly  wide  range 
of  lift  coefficients  and  thickness-to-chord  ratios,  but  no  data 
were  available  giving  the  basic  characteristics  of  the  section 
equipped  with  a  control  flap  for  hydrofoil  use.  Such  data  are 
desirable  for  hydrofoil  ship  design  and  performance  predictions 
and  can  provide  a  valuable  check  on  the  accuracy  of  prediction 
techniques  being  developed  for  the  analysis  of  flow  over 
flapped  hydrofoils.  They  can  also  be  expected  to  provide  a 
standard  for  comparison  with  new  hydrofoil  sections  developed 
to  take  advantage  of  improved  analytical  techniques. 

To  meet  the  requirement  for  basic  data,  a  program  of 
tests  was  undertaken  on  a  NACA  16-309  section  with  an  a  =  1.0 
mean  camber  line  and  with  a  25%  flap-chord  ratio,  simple,  sealed 
flap.  This  particular  section  was  selected  because  it  is  fitted 
on  USS  "HIGHPOINT"  (PCH-1),  a  128  ton  auw  hydrofoil  ship  used 
extensively  by  DTNSRDC  for  hydrofoil  research.  It  consequently 
presents  the  best  opportunity  for  comparison  with  three- 
dimensional  model  tests  and  with  full  scale  operation. 

The  section  test  program  included  both  water  tunnel 
and  wind  tunnel  tests.  The  primary  objective  of  water 
tunnel  tests,  carried  out  at  the  California  Institute  of  Tech¬ 
nology  (CIT),  was  to  determine  the  cavitation  characteristics 
for  a  range  of  flap  angles  and  angles  of  attack,  and  their 
effects  on  the  lift,  drag  and  pitching  moment  characteristics. 
The  wind  tunnel  tests  were  made  to  obtain  more  accurate  force 
and  moment  data  and  also  to  obtain  the  pressure  distribution 
characteristics,  since  these  are  obtained  much  more  readily 
in  the  wind  tunnel. 
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The  purpose  of  this  report  is  to  present  and  analyze 
the  data  from  the  water  tunnel  tests.  The  model,  experimental 
facility,  instrumentation  and  test  procedure  are  described  in 
some  detail.  A  fairly  extensive  analysis  of  the  experimental 
data  is  then  presented,  from  which  a  number  of  significant  con¬ 
clusions  are  drawn.  For  completeness,  the  test  data  plots  are 
given  in  full  in  Appendix  A. 

A  consideration  of  fundamental  importance  in  any  model 
test  program  is  the  influence  of  scale  effects.  Recent  work2’3 
indicates  that  cavitation  inception  is  not  simply  dependent  on 
the  local  static  pressure.  It  can  originate  as  a  result  of 
fairly  large  pressure  fluctuations  associated  with  transition 
from  laminar  to  turbulent  flow.  The  location  of  the  transition 
point  along  the  length  of  the  chord  is  Reynolds  number  dependent 
and  thus,  the  inception  of  cavitation  can  also  be  Reynolds  num¬ 
ber  dependent.  The  CIT  tests  were  made  mostly  at  a  tunnel 
water  speed  of  50  ft/sec,  giving  a  chord  Reynolds  number  of 
2.49  x  106.  The  highest  Reynolds  number  possible  was  3.98  x 
10°  which,  although  high  by  model  test  standards,  is  still 
significantly  smaller  than  the  representative  value  for  USS 
HIGHPOINT  of  2.5  x  107.  Tests  were  therefore  made  to  investi¬ 
gate  the  sensitivity  of  cavitation  occurrence  to  Reynolds  number 
over  the  available  range.  Reynolds  number  sensitivity  tests 
were  also  made  on  the  force  and  moment  data. 

The  wind  tunnel  test  series'*,  which  was  run  on  a 
similar  model  of  30  inches  chord  at  Reynolds  numbers  up  to 
4  x  10*’,  provides  further  data  on  scale  effects.  Use  is  made 
of  these  data  in  discussing  scale  effects  arising  from  different 
boundary  layer  flow  characteristics. 


2,  MODEL  DESCRIPTION 


The  six  inch  chord,  six  inch  span  water  tunnel  model 
was  made  from  17-4PH  steel  hardened  to  the  H900  condition.  It 
is  shown  in  Figure  1,  complete  with  tunnel  wall  fairing  disc 
and  ready  for  installation.  The  flap  could  be  positioned  at 
discrete  angle  settings  of:  -10°,  -5°,  0°,  2.5°,  5°,  7.5°,  10°, 
15°  and  20°.  The  setting  plate  and  the  hole  alignments  used 
to  obtain  the  required  flap  angles  are  illustrated  in  Figure  2 
which  shows  the  model  with  fairing  disc  removed.  A  separate 
securing  plate  was  required  at  the  free  e  d  of  the  model  for 
each  flap  angle  setting. 

The  section  coordinates  were  derived  from  the  standard 
data  given  by  Abbott  and  Von  Doenhoff5,  using  a  cubic  spline  fit 


info  a  leading  edge  radius  of  0.024  inch.  The  coordinates  are 
given  in  TabLe  I.  Profile  measurements  were  made  on  the  com¬ 
pleted  model  by  CIT,  along  the  axes  shown  in  Figure  3.  Measure¬ 
ments  along  axes  A,  B  and  C  indicate  virtually  no  twist  along 
the  span.  Deviations  from  the  specified  section  profiles  along 
D,  E  and  F  are  illustrated  in  Figure  4.  This  shows  the  error- 
to-chord  ratio  to  be  generally  less  than  0.0005.  In  the  vicinity 
of  the  hinge-line,  there  is  a  deviation  in  the  surface  contour 
which  is  higher  and  steeper  than  desirable.  There  is  a  second 
deviation  at  about  5%  chord,  not  as  large  in  amplitude  but 
equally  steep.  The  model  finish  was  generally  fine  and  polished. 
Castings  of  the  leading  edge,  sectioned  at  mid-span,  were  in¬ 
spected  on  an  optical  comparator  at  a  magnification  of  20X  and 
revealed  no  roughness  or  discontinuity  in  the  curvature.  It  is 
difficult  to  obtain  very  high  profile  accuracies  in  the  small 
size  of  model  required  for  water  tunnel  testing.  These  profile 
errors  are  considered  to  be  within  the  normal  range  of  expect¬ 
ation  for  a  relatively  complicated  model  incorporating  a  hinged 
flap. 


3.  FACILITY  AND  INSTRUMENTATION 


The  High  Speed  Water  Tunnel  in  the  Graduate  Aero¬ 
nautical  Laboratories  of  the  California  Institute  of  Technology 
(GALCIT)'3  was  used  in  the  performance  of  this  experiment.  For 
these  tests,  the  tunnel  was  equipped  with  the  two-dimensional 
working  section  shown  in  Figure  5.  Test  section  pressure  and 
water  velocity  could  be  controlled  independently  over  a  wide 
range  and  gases  released  by  cavitation  on  the  mode  1  were  re¬ 
absorbed  before  the  flow  reentered  the  working  section. 

The  model  was  viewed  through  a  four- inch  thick  " P L e x  i  - 
glas"  window,  the  inside  surface  of  which  was  parallel  to  the 
theoretical  centerline  of  the  tunnel.  The  opposite  wall  was 
adjustable  and  rigidly  attached  to  the  tunnel  at  the  upstream 
end.  The  position  of  this  wall  was  adjusted  to  compensate  for 
the  effects  of  longitudinal  pressure  gradient  due  to  wall 
boundary  layer  growth.  The  model  fairing  plate  was  fitted  flush 
with  the  wall  and  great  care  was  taken  to  avoid  spurious  flow 
disturbances  due  to  discontinuities. 

The  hydrofoil  model  was  supported  as  a  cantilever  from 
a  three-component  strain  gage  balance  equipped  with  interchange¬ 
able  load  cells.  The  highest  ranges  available  were  used  for 
these  tests.  These  were:  lift  +  800  lbs,  drag  +  200  lbs,  pitch¬ 
ing  moment  +  1200  inch-lbs.  Each  of  the  three  output  signals 
was  amplified,  digitised  and  simultaneously  time-averaged. 
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Signals  were  displayed  on  panel  meters  and  recorded  on  punched 
cards. 


Working  section  velocity  was  determined  by  measuring 
the  static  pressure  difference  between  the  settling  chamber 
upstream  of  the  nozzle  and  the  entrance  to  the  working  section, 
using  a  mercury-water  manometer.  Manometer  readings  were  con¬ 
verted  to  velocity  through  an  experimentally-determined  cali¬ 
bration  factor,  based  on  measurements  made  with  a  Prandtl  type 
pitot-static  tube.  Working  section  pressure  was  measured  by 
means  of  a  mercury  barometer  connected  to  the  piezometer  row 
at  the  entrance  to  the  working  section  and  referenced  to  the 
horizontal  centerline  of  the  working  section. 

The  axial  turbulence  level,  measured  on  the  centerline 
at  the  balance  station  was  less  than  0.25%  and  was  estimated  by 
CIT*'  to  be  0.2%.  The  dissolved  gas  content  was  measured  at  be¬ 
tween  12  and  14  parts  per  million. 


4.  PROCEDURE 


The  model  was  mounted  wit^  its  chord  line  parallel  to 
the  centerline  of  the  working  section  for  an  angle  of  attack  of 
zero  degrees.  The  working  section  wall  was  set  to  give  a  nearly 
zero  horizontal  pressure  gradient  over  the  range  of  water  vel¬ 
ocities. 


At  the  start  of  the  program,  the  conventional  "Plexi¬ 
glas"  viewing  window  was  used,  with  the  gap  between  it  and  the 
free  end  of  the  model  set  to  0.0)0  inch.  This  gap  was  required 
to  prevent  Interference  as  the  window  deflects  under  differential 
pressure  when  tile  working  section  pressure  was  reduced.  For  the 
later  part  of  the  program,  the  "Plexiglas"  sidewall  was  replaced 
by  a  stronger,  aluminum  plate,  1.75  inches  thick,  with  a  small, 
7.5  inch  diameter  viewing  port.  The  tip  gap  was  reduced  to 
0.020  inch  for  this  installation.  The  type  of  window  used  is 
listed  in  Table  II  for  the  given  test  conditions.  The  aluminum 
plate  was  used  for  all  Reynolds  number  dependence  tests,  where 
positive  tank  pressures  were  involved  and,  as  a  matter  of  con¬ 
venience,  for  all  tare  runs.  Measurements  of  window  deflection 
as  a  function  of  working  section  pressure  are  presented  in  Fig¬ 
ure  6 . 


Since  the  model  had  to  be  removed  from  the  tunnel 
each  time  that  the  flap  angle  was  adjusted,  test  series  were 
run  with  fixed  flap  angle  and  with  angle  of  attack,  Reynolds 
or  cavitation  numbers  as  variables.  For  the  Reynolds  number 
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surveys,  runs  were  made  at  fixed  angle  of  attack  and  flap  angle 
at  speeds  of  25,  40,  50,  70  and  80  ft/sec,  giving  chord  Reynolds 

numbers  of  1.25,  1.99,  2.49,  3.49  and  3.98  x  106  respectively. 

Tunnel  pressures  were  increased  or  decreased  as  required  to 
obtain  cavitation-free  running  and  to  measure  the  cavitation 
numbers  for  inception  and  desinence  at  the  leading  edge  and 
hinge-line.  The  corresponding  balance  measurements  were  then 
recorded . 


For  the  cavitation  survey,  runs  were  made  at  a  con¬ 
stant  water  speed  of  50  ft/sec,  a  compromise  speed  which  al¬ 
lowed  reasonable  ranges  of  flap  angle  and  angle  of  attack,  while 
generally  keeping  within  balance  loading  restrictions.  In  each 
case,  force  data  were  measured  at  cavitation  numbers  of  2.04, 
1.26,  0.90,  0.53,  0.33  and  0.22.  Data  were  also  taken  for  the 
cavitation  numbers  associated  with  incipient  and  desinent  cavit¬ 
ation. 


When  the  cavity  extended  past  mid-chord,  large  unsteady 
forces  developed,  with  peaks  observed  by  oscilloscope  to  be  two 
to  four  times  steady  state.  These  overloaded  the  balance  and 
made  it  impossible  to  make  measurements  over  the  wide  ranges  of 
flap  angle  and  angle  of  attack  originally  planned.  An  additional 
restriction  occurred  at  high  flap  angles,  when  the  working  section 
pressures  required  to  suppress  cavitation  exceeded  the  maximum 
a  1 lo wab 1 e . 


5,  CORRECTIONS  AND  ACCURACY 


The  force  balance,  manometer,  physical  constants  and 
other  data  were  recorded  on  punched  cards  and  reduced  using 
digital  computer  processing. 

The  balance  zero  values  were  dependent  on  the  model 
attitude,  making  corrections  necessary  for  these  "gravity  tare" 
values  for  each  run.  The  measured  data  were  also  corrected  for 
the  effect  of  balance  interactions  bv  the  application  of  an  in¬ 
verse  calibration  matrix.  Finally,  the  data  were  corrected  for 
the  influence  of  the  forces  acting  on  the  fairing  disc,  the 
"balance  tares".  These  were  determined  for  each  model  test  con¬ 
dition  using  the  following  procedure.  The  model  was  removed 
from  the  force  balance,  inverted  and  installed  in  a  support 
fixture  on  the  opposite  wall.  A  dummy  fairing  disc  was  then 
mounted  on  the  force  balance  and  the  tip  of  the  model  brought 
to  within  0.002  inch  of  the  disc.  This  clearance  was  maintained 
by  mechanically  adjusting  the  hydrofoil  model  to  compensate  for 
changes  due  to  working  section  pressure.  The  disc  forces  could 
then  be  measured  for  each  test  condition  in  flow  conditions 
similar  in  many  respects  to  those  applying  in  practice. 


5 


I 


A 


—  *, 


nr.1  '•  myi  jt  y. .  >  j 


-=■ 


These  balance  Care  values  were  reasonable  for  lift 
and  moment,  representing  at  their  worst,  less  than  1  %  and  3% 
of  the  maximum  recorded  values.  This  was  not  the  case  for  the 
drag  measurement.  Maximum  drag  tare  values  were  about  25%  of 
maximum  recorded  drag  and  appear  to  have  varied  systematically 
with  tunnel  pressure.  This  is  illustrated  in  Figure  7  which 
shows  the  drag  tare  values  for  the  0°a,  0°5  and  0°a»  5°6  cases 
plotted  against  tunnel  pressure  for  various  speeds.  The  tare 
measurement  arrangement  was  a  close  approximation  to  and  not 
an  exact  duplication  of  measurement  conditions,  and  it  is  difficult 
to  say  what  influences  were  present  in  practice.  The  fact  that 
drag  tare  values  could  be  large  and  variable  is  considered  to 
prejudice  significantly  the  reliability  of  the  drag  results. 

No  tunnel  boundary  corrections  were  applied  and  there 
were  no  corrections  for  model  deflections  under  load. 


Manometer  indications  of  velocity  in  the  working 
section  exhibited  a  fluctuation  which  varied  from  approximately 
+  0.5%  at  10  ft/sec  to  +  1%  at  60  ft/sec.  The  period  was  long 
and  there  were  frequent  occasions  when  the  reading  was  virtually 
steady.  Data  were  recorded  during  these  occasions  and  as  a 
result,  velocity  data  presented  here  are  believed  accurate  to 
within  +  0.5%.  Working  section  pressure  indicators  exhibited 
similar,  though  smaller,  fluctuations  and  these  measurements 
are  also  believed  accurate  to  within  4-  0.5%. 


CIT  estimates  for  other  total  measurement  errors  are 

as  f  ollows : 

angle  of  attack:  +  0.1  deg 


data  reduction : 
lift  force: 


+  0.1%  of  indicated  value 

+  2.1  lbs,  giving 

Cn  =  +  .0035  at  50  ft/sec 


pitching  moment:  +  3.1  in-lbs,  giving 


C  =  +  .0051  at  50  ft/sec 
m  — 


6.  DISCUSSION  OF  RESULTS 


6.1  General  Note  on  Boundary  Layer  Flow  Conditions 

and  Effects 


The  NACA  16  Series  section  is  designed  to  have  an 
extensive  region  of  laminar  flow  for  a  range  of  lift  coefficients 
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around  the  design  value,  giving  comparatively  low  drag  under 
normal  steady  state  conditions.  The  location  of  the  point  of 
transition  to  turbulent  flow  moves  forward  suddenly  when 

is  either  decreased  or  increased  beyond  the  laminar  range, 
leading  to  a  sharp  increase  in  drag.  Flow  is  then  said  to  have 
become  turbulent.  The  range  of  values  for  which  laminar 

flow  conditions  are  realised  becomes  narrower  with  increasing 
Reynolds  number. 

The  companion  series  of  wind  tunnel  tests4  shed  con¬ 
siderable  light  on  the  boundary  layer  conditions  applying  to  the 
wind  and  water  tunnel  models.  Laminar  flow  conditions  were 
obtained  as  expected  in  the  wind  tunnel  but  there  was  an  ad¬ 
ditional  complication  because  transition  to  turbulent  flow  was 
accompanied  by  flow  separation  over  the  trailing  edge  of  the 
flap.  This  in  turn  led  to  sharp  decreases  in  the  effective 
camber  of  the  section,  decreasing  the  section  lift.  The  largest 
decreases  occurred  for  higher  flap  angles  at  the  positive  limit 
of  the  laminar  Cr,  range  but  similar,  smaller  changes  occurred 

at  the  negative  limit. 


Wind  tunnel  lift  coefficients  arc  shown  plotted  against 
angle  of  attack  in  Figure  8,  both  with  and  without  roughness 
strips.  Application  of  roughness  strips  to  the  upper  and  lower 
surfaces  of  the  section  at  the  5%  chord  points  fixed  the  points 
of  transition  over  the  whole  range  and  established  the  turbu¬ 
lent  flow  condition  artificially.  This  changed  the  nature  of  the 
force  and  moment  characteristics  because  trailing  edge  separation 
then  took  place  comparatively  smoothly  with  increasing  and  de¬ 
creasing  flap  angle  and  angle  of  attack.  The  lift  values  were 
reduced  over  the  formerly  unseparated  range  but  the  sharp  dis¬ 
continuities  associated  with  delayed  transition  to  turbulent 
flow  were  eliminated. 


The  lift  data  from  the  water  tunnel  are  also  shown  in 
Figure  8.  These  have  values  which  lie  between  the  transition- 
fixed  and  t r an s i t ion- f r e e  wind  tunnel  data,  over  the  formerly 
unseparated  part  of  the  wind  tunnel  range  of  angles  of  attack. 
The  indications  are  that  for  the  water  tunnel  model,  separation 
spread  smoothly  over  the  trailing  edge  with  increasing  angle  of 
attack  in  a  manner  somewhat  similar  to  the  t r ans i t i on- f ixed  wind 
tunnel  model.  The  reason  for  this  is  not  clear  since  the  free 
stream  turbulence  levels  are  thought  to  be  similar  in  both 
facilities  and  chord  Reynolds  number  was  lower  in  the  water 
tunnel.  However  the  surface  roughness  Reynolds  number  may  well 
have  been  higher  for  the  water  tunnel  model.  This  and  the 
relatively  large  section  profile  errors  could  have  been  factors 
in  promoting  more  rapid  transition  and  spread  of  trailing  edge 
separation. 
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6.2  Lift  Characteristics 


Lift  coefficient  is  plotted  against  angle  of  attack 
in  Figure  9  for  the  full  range  of  flap  angles  and  for  a  cavit¬ 
ation  number  of  2.04,  representing  the  largely  uncavitated  con¬ 
dition.  The  curves  are  linear  for  small  ranges  of  angle  of 
attack,  and  flap  angle,  but  show  a  trend  to  increasing  non¬ 
linearity  as  angle  of  attack  and  flap  angle  approach  the  test 
range  limits.  The  zero-lift  angle  of  -2.1°  for  0°  flap  angle 
is  less  negative  than  that  obtained  in  both  the  available  sets 
of  wind  tunnel  data1*1*  which  gave  a  value  of  -2.6°  for  the 
t r ans i t ion- f r ee  test  condition.  The  non-linearity  of  the  data 
and  less  negative  zero-lift  angle  are  both  attributed  to  flow 
separation  over  the  trailing  edge  of  the  flap,  as  discussed 
earlier  in  Section  6.1. 

Lift  coefficient  is  shown  plotted  against  flap  angle 
with  angle  of  attack  as  a  parameter  in  Figure  10,  for  the  same 
cavitation  number  of  2.04.  The  curves  depart  substantially 
from  linear  for  both  positive  and  negative  flap  angle,  again 
apparently  as  a  result  of  flow  separation  on  the  upper  or  lower 
surface  of  the  flap. 

The  typical  effect  of  cavitation  number  on  lift  co¬ 
efficient  at  0°  angle  of  attack  is  shown  in  Figure  11  for  the 
full  range  of  flap  angles.  The  cavitation  inception  limit  lines 
on  the  diagram  show  that  leading  edge  cavitation  could  be  fol¬ 
lowed  by  an  increase  in  lift  as  cavitation  number  decreased 
(see  also  Figures  A13  to  A19  in  Appendix  A).  Presumably  the 
effect  of  a  small  amount  of  leading  edge  cavitation  was  to 
increase  the  effective  section  camber.  As  the  cavitation  number 
of  the  flow  was  decreased,  cavitation  spread  chordwise  over 
the  section  until  it  reached  the  mid-back  or  hinge-line  when 
the  forces  tended  to  become  very  unsteady.  Marked  decrease  of 
lift  coefficient  was  normally  associated  with  extension  of 
hinge-line  or  leading  edge  cavitation  beyond  the  trailing 
edge,  when  the  effective  camber  would  be  much  decreased. 

The  travelling,  growing  bubble  type  of  cavitation 
was  the  type  most  often  observed  but  at  high  angles  of  attack 
or  flap  angle,  sheet  cavitation  could  occur  in  a  band  at  the 
leading  edges  of  both  the  foil  and  flap. 

The  foil  and  flap  lift-curve  slopes  are  plotted 
against  flap  angle  in  Figure  12.  These  are  defined  as  the 
rates  of  change  of  lift  coefficient  with  angle  of  attack  and 
flap  angle  respectively  and  are  shown  for  the  uncavitated  con¬ 
dition  at  0°  angle  of  attack.  The  foil  lift-curve  slope  has 
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a  value  of  0.096  per  degree  for  a  small  range  of  flap  angles 
close  to  0°  (about  4%  less  than  the  corresponding  wind  tunnel 
value)  but  falls  off  sharply  for  flap  angles  greater  than  5°. 

The  flap  lift-curve  slope  has  a  maximum  value  of  about  0.06 
for  a  very  small  range  of  flap  angles  close  to  zero. 

The  corresponding  flap  effectiveness  ratios,  defined 
as  the  ratio  of  flap  lift-curve  slope  to  foil  lift-curve  slope, 
are  plotted  in  Figure  13  for  a  constant  lift-curve  slope  of 
0.096  per  degree.  The  theoretical  value  can  be  approximated 
for  the  unseparated  flow  condition  by7 

^  =  -/f~=  0.64 

C*a  U 

This  value  is  achieved,  but  only  for  a  very  small  range  of  flap 
angles  close  to  zero. 

6 . 3  Drag  Characteristics 

Drag  coefficient  is  plotted  against  lift  coefficient 
in  Figure  14  for  a  range  of  flap  angles  at  a  cavitation  number 
of  2.04.  Data  from  the  companion  wind  tunnel  tests,  included 
for  comparison,  show  the  water  tunnel  data  to  be  higher  in 
value  except  for  a  small  range  of  values  close  to  the  design 

point.  The  wind  tunnel  values  were  obtained  for  a  large,  ac¬ 
curate  model,  using  the  comparatively  precise  wake  survey  method. 
Despite  the  care  taken  with  the  water  tunnel  force  and  moment 
measurements,  the  balance  and  the  model  mounting  arrangements 
seem  likely  to  have  introduced  errors,  as  indicated  by  the  com¬ 
paratively  large  drag  tare  values  mentioned  in  Section  5.  Con¬ 
sequently,  with  the  exception  of  Figures  14  and  19  no  drag  data 
are  included  here. 

6 . 4  Moment  Characteristics 

The  quarter-chord  pitching  moment,  C  ,  is  plotted  as 

m 

a  function  of  lift  coefficient  in  Figure  15  for  the  range  of 
flap  angles  at  a  cavitation  number  of  2.04.  C^  can  be  expected 

to  become  more  negative  as  flap  angle  increases  owing  to  the 
increase  in  effective  camber.  The  rate  of  change  can  be  estim¬ 
ated  from  the  thin  aerofoil  theory  of  Glauert8  as 

C  =  -2/f (1  -  f ) 5  =  -0.011  per  degree, 
m  o 

Figure  15  shows  that  a  C  change  of  approximately  0.055  results 

m 

from  a  5°  flap  angle  change  at  values  of  C^  close  to  0.2.  How- 
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ever,  there  is  also  a  marked  change  of  C  with  C  „  .  This  is  pre- 

m  Jc 

sumably  due  to  increasing  separation  and  straightening  of  the 
flow.  Because  of  this,  a  22%  chord  location  for  the  aerodynamic 
center  appears  to  be  a  better  compromise  than  the  standard  value 
of  25%  for  an  unflapped  section. 

Figure  16  shows  the  effect  of  flow  cavitation  number 
on  pitching  moment  coefficient  for  0°  angle  of  attack.  The 
data  are  generally  consistent  with  the  lift  data,  increasing 
sharply  as  cavitation  number  decreases  to  very  low  numbers,  but 
not  necessarily  changing  immediately  with  the  onset  of  cavitation. 

6.5  Reynolds  Number  Dependence  of  Force  and  Moment 
Coef  f icients 


The  variation  of  lift  with  Reynolds  number  is  shown 
In  Figure  17  for  the  four  foil  and  flap  angle  combinations  given 
by  a  =  0°  and  2°  and  6=0°  and  5°.  The  lift  values  for  the 
08a,  5°6  case  show  a  comparatively  sharp  decrease  with  increas¬ 
ing  Reynolds  number.  This  is  seen  from  Figure  8  to  be  a  likely 
result  of  increased  flow  separation  at  the  flap  trailing  edge. 

There  is  little  change  at  the  other  angle  settings. 
Figure  8  shows  that  any  change  due  to  increasing  separation  at 
the  0°  flap  angle  settings  would  be  comparatively  small.  At  2°a, 
5°6,  flow  separation  is  apparently  well  established  even  at  the 
lowest  Reynolds  number  setting. 

The  corresponding  quarter  chord  pitching  moment  data 
of  Figure  18  are  generally  consistent  with  this  analysis.  There 
is  a  small  increase  In  the  0°  flap  angle  data  values  with  Rey¬ 
nolds  number,  while  the  <)°ii,  5°,S  values  approach  those  for  the 

2°u,  5°iS  setting.  This  ;s  to  be  expected  since  quarter  chord 
pitching  moment  coefficients  ere  dependent  on  effective  section 
camber  and  will  be  similar  foi  the  same  flap  angle,  when  the 
boundary  layer  conditions  are  similar. 

The  variation  of  drag  coefficient  with  Reynolds  number 

is  shown  in  Figure  19  for  the  0°a,  0°6  case,  the  only  angle 

setting  combination  for  which  the  drag  data  may  be  reasonably 

representative.  The  C,  values  increase  with  Re  for  Re  <  2.49 

a 

x  106  and  then  tend  to  follow  the  well-established  trend  line7 
for  smooth  surfaces  in  turbulent  flow,  indicating  that  at  this 
angle  setting,  transition  to  turbulent  flow  occurs  at  about 
Re  »  2.49  x  10  6 . 
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6 . 6  Reynolds  Number  Dependence  of  Cavitation 


Tests  were  made  for  cavitation  dependence  on  Reynolds 
number  at  the  same  four  angle  settings,  using  the  procedure 
described  in  Section  4.  The  flow  cavitation  numbers  at  which 
inception  and  desinence  occurred  are  plotted  in  Figure  20  for 
the  0°a,  0° 5  setting.  For  this  case,  cavitation  took  place  virtu¬ 
ally  simultaneously  on  the  upper  surface  at  the  flap  hinge-line 
and  at  section  mid-back,  with  little  hysteresis  ard  virtually 
no  variation  with  Reynolds  number. 

For  the  a  =  0°,  5=5°  case,  shown  in  Figure  21,  hinge¬ 

line  cavitation  occurred  first  as  cavitation  number  was  decreased, 
followed  for  Re  <  2.49  x  106  by  mid-back  cavitation.  For  Re  > 

2.49  x  106,  however,  although  hinge-line  cavitation  still 
occurred  first,  it  was  followed  by  leading  edge  cavitation. 

For  this  particular  setting,  wind  tunnel  data  show  the  pressure 
distribution  to  be  very  flat  over  the  first  half  of  the  chord, 
so  that  little  change  would  be  required  in  the  flow  to  alter  the 
cavitation  characteristics.  In  fact,  the  lift  and  moment  co¬ 
efficient  values  undergo  fairly  large  changes  with  Reynolds  number, 
as  noted  in  Section  6.5  above,  indicating  a  significant  change  to 
the  flow  pattern.  In  these  circumstances,  the  switch  from  mid- 
back  to  leading  edge  cavitation  is  not  surprising. 

Cavitation  took  place  only  at  the  upper  surface  lead¬ 
ing  edge  for  the  a  =  2°,  5=0°  setting,  as  shown  in  Figure  22. 

Some  of  the  data  exhibit  considerable  hysteresis  and  there  are 
large  differences  between  the  data  points  obtained  during  the 
Reynolds  number  survey  and  those  obtained  during  the  general 
survey  at  Re  =  2.49  x  106.  There  are  too  few  data  points  to 
establish  the  trends  with  certainty  but  the  cavitation  inception 
number  appears  to  switch  in  value  over  the  Reynolds  number  test 
range,  with  some  uncertainty  at  intermediate  Re  values. 

Cavitation  also  occurred  only  at  the  upper  surface 
leading  edge  for  the  a  =  2°,  6*5°  setting  and  Figure  23  in¬ 
dicates  a  strong  Reynolds  number  dependency.  For  this  setting, 
cavitation  spread  chordwise  with  decreasing  Reynolds  number 
until  it  reached  the  mid-back  when  it  spontaneously  oscillated 
between  there  and  the  trailing  edge  of  the  flap.  The  resulting 
violent  force  fluctuations  prevented  testing  over  the  full 
Reynolds  number  range. 

These  tests  for  Reynolds  number  dependence  of  cavit¬ 
ation  indicate  that  mid-back  and  hinge-line  cavitation  are  in¬ 
dependent  of  Reynolds  number  over  the  test  range.  Leading 
edge  cavitation  occurrence  does  seem  to  be  Reynolds  number 
dependent  but  there  are  insufficient  data  to  assess  the  extent 
to  which  this  dependence  is  systematic. 
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The  limit  lines  for  cavitation  inception  are  given 
in  Figure  24  for  all  flap  angles  tested.  The  hinge-line  (or 
occasionally  mid-back)  limits,  represented  by  the  lower  lines 
of  the  bucket  curves,  increase  with  flap  angle  as  would  be  ex¬ 
pected.  The  Reynolds  number  dependency  tests  described  in 
Section  6.6  indicate  that  these  lines  will  be  relatively  in¬ 
dependent  of  Reynolds  number. 


The  leading  edge  cavitation  limits  (represented  by 
the  right  hand  lines)  on  the  other  hand,  will  be  affected  by 
Reynolds  number.  They  seem  likely  to  move  to  the  left  with 
increasing  Reynolds  number,  narrowing  the  range  of  uncavitated 
available  at  the  higher  full  scale  Reynolds  numbers. 


7.  CONCLUSIONS 

Tests  on  a  16-309  section  equipped  with  a  25%  flap 
chord  ratio,  simple,  sealed  flap  show  the  characteristics  to 
be  influenced  greatly  by  turbulent  flow  separation  over  the 
flap  trailing  edge.  This  greatly  decreases  the  lift  available 
over  much  of  the  normal  range  of  flap  angle  and  angle  of  attack. 

Values  of  foil  lift-curve  slope  and  flap  effective¬ 
ness  are  close  to  theoretical  at  0°  flap  angle  but  decline 
fairly  rapidly  as  flap  angle  increases  or  decreases. 

The  drag  coefficient  data  are  unduly  high  and  are  con¬ 
sidered  suspect  except  for  a  narrow  range  of  lift  coefficients 
between  about  zero  and  0.25.  The  reason  for  this  is  not  known 
but  may  have  been  due  to  difficulties  in  accounting  for  balance 
tares . 

The  pitching  moment  data  indicate  the  aerodynamic  center 
to  be  at  about  0.22  chord.  For  flap  angles  greater  than  5°, 
considerable  increase  in  the  quarter  chord  pitching  moment  co¬ 
efficient  occurs  with  increasing  C^. 

Comparison  with  companion  wind  tunnel  test  data,  to¬ 
gether  with  Reynolds  number  sensitivity  tests  indicate  that 
transition  to  turbulent  flow  occurs  fairly  readily  in  the  water 
tunnel.  Turbulent  flow  is  well  established  on  the  model  at  the 
main  test  Reynolds  number  of  2.49  x  106  but  the  indications  are 
that  some  forward  movement  of  the  transition  point  continues 
as  Reynolds  number  is  increased. 

Despite  well-established  turbulent  flow,  leading 
edge  cavitation  occurrence  is  sensitive  to  Reynolds  number, 
making  the  data  useful  only  as  a  general  guide  to  full-scale 
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performance.  Mid-back  and  hinge-line  cavitation  proved  to  be 
insensitive  to  Reynolds  number  over  the  test  range. 

For  the  higher  angle  of  attack  cases,  where  leading 
edge  cavitation  occurs  well  ahead  of  mid-back  or  hinge-line 
cavitation  as  cavitation  number  decreases,  leading  edge  cavit¬ 
ation  does  not  lead  to  immediate  loss  of  lift  or  increase  of 
drag.  Marked  force  and  moment  changes  are  more  normally 
associated  with  the  extension  of  the  mid-back  or  hinge-line 
cavity  beyond  the  trailing  edge. 
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